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Abstract 

The hyperfine interaction of phosphorus donors in fully strained Si thin films grown on virtual 
Sii-^Gea; substrates with x < 0.3 is determined via electrically detected magnetic resonance. For 
highly strained epilayers, hyperfine interactions as low as 0.8 mT are observed, significantly below 
the limit predicted by valley repopulation. Within a Green's function approach, density functional 
theory (DFT) shows that the additional reduction is caused by the volume increase of the unit cell 
and a local relaxation of the Si ligands of the P donor. 
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At present, several approaches for solid-state based quantum computing hardware are 
actively pursued. Thepossible integration with existing microelectronics and the long de- 
coherence times are particular advantages of concepts [4j, |5j, |g, |7[ using the nuclear 
or electronic spins of phosphorus donors in group IV semiconductors as qubits. These con- 
cepts require gate-controlled exchange coupling between neighbouring donors. However, to 
control the exchange coupling in semiconductors with an indirect band structure such as 



Si, the donor atoms have to be positioned with atomic precision 



due to Kohn-Luttinger 



oscillations of the donor wavefunction 0, IJQ]. Under uniaxial compressive strain in [001] 
direction, two conduction band (CB) minima are lowered in energy, which leads to a sup- 
pression of the oscillatory behaviour of the wavefunction for donors located in the (001) 
lattice plane |9(. The strain will also affect the wavefunction at the position of the donor 
atom, which can be observed directly via the hyperfme (hf) interaction between the donor 



electron and its nucleus 
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. In this letter, we experimentally and theoretically study the 



hf interaction of phosphorus donors in silicon as a function of uniaxial compressive strain in 
thin layers of Si on virtual SiGe substrates, extending the regime investigated by Wilson and 
Feher |ll| by a factor of 20 to higher strains. We find that the reduction of the hf interaction 
significantly exceeds the limit predicted so far 

fly 

.lll|. By ab-initio DFT calculations using 



a Green's function approach 



12|, we are able to show that this reduction is caused by the 



strain-induced lattice distortion, and a relaxation of the Si lattice next to the P atoms. As 
the nuclear exchange mediated by the overlap of the donor wavefunction scales with the 
square of the hf interaction, these results have a direct impact on the rate at with two-qubit 
operations can be performed jj]. 

The samples under study were prepared by chemical vapor deposition (CVD) on 30 ficm 
B-doped Si (001) substrates jl^ . Fully strained thin P-doped Si epilayers with [P] ~ 
1 x 10 17 cm -3 were grown lattice matched on virtual relaxed Sii^Ge^-substrates optimized 
for low dislocation densities. These substrates consist of an intrinsic 0.3 //m Si buffer, a 
variable sequence of Sii-^Ge^, buffer layers with stepwise increasing y, and the actual 2 /iin 
thick virtual Sii^Ge^,. substrate with Ge-contents x = 0.07,0.15,0.20,0.25, and 0.30. High- 
resolution X-ray diffraction (HRXRD) on the 004 and 224 reflexes was used to determine 
the (in- and out-of-plane) lattice constants of the SiGe layers from which the Ge content and 
degree of relaxation were calculated using a parabolic dependency of x on the relaxed lattice 
constants ^| and a linear interpolation of the elastic stiffness constants between Si and Ge 
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0. 



[15(. The Sii-aGe^ layer determines the strain of the Si epilayer: The larger lattice constant 
of SiGe alloys compared to Si leads to biaxial tensile strain, accompanied by a compensating 
uniaxial compressive strain in growth direction as shown by the inset in Fig. ^ All Si:P 
epilayers have a thickness of 15 nm, well below the critical thickness of Si on Sio.7Geo.3 (l6|. 

Figure n shows the corresponding HRXRD reciprocal space-map (RSM) of the 224 reflex 
for the Sio.7Geo.3 sample. T he p eak arising from the Si wafer is set to match the lattice 
constant of c-Si (5.4310 A) |l7[. Extending to lower reciprocal lattice units q\\=y/2\/ a\\ 
(in plane) and q^=\fA\/ f a± (out of plane), a set of diffraction peaks arising from the SiGe 
buffer layers is observed, followed by the peak of the virtual substrate. The slight difference 
between the observed values and those expected for fully relaxed Sii_ x Ge x alloys (black solid 
line, ^j]) indicates a degree of relaxation of 96.4% of the virtual substrate. Additionally, 
Fig.dshows at g ± =0.5715 (a ± =5.3912 A) and g||=0.39686 (a||=5.4898 A) the diffraction by 
the strained silicon epilayer. The perfect agreement of qu of the Si epilayer and the virtual 
substrate reflects the fully strained pseudomorphic growth. It is important to note that 
the compression in growth direction is in accordance with linear elasticity theory. With 
Cxi = 165.6 GPa and Cu = 63.9 GPa as the stiffness constants of Si 15| . the relation 
e_i_=— 2C12/C11 • 6|| between in-plane and out-of-plane strains e\\ = (a^-a Si )/a Si and e_|_=(a_i_- 
a sO/ a Si predicts a value of a_i_=5.3856 A on the Si .7Ge .3 substrate, in reasonably good 
agreement with a± obtained from Fig. ^ taking into account the width of the thin Si layer 
caused by the small thickness. Hence, the strain values for the different strained Si samples 
are obtained from the relaxed in-plane lattice constants of the corresponding SiGe buffers 
using linear elasticity theory. 

To observe the the small amount of P donors in the thin strained layer with high sensitiv- 
ity, we use electrically detected magnetic resonance (EDMR), which monitors spin resonance 
via the influence of spin selection rules on charge transport processes j^, 19, 20]. The EMDR 
experiments were performed in a dielectric ring microwave resonator using a HP83640A mi- 
crowave source and measuring the spin-dependent photoconductivity at T = 5 K in a liquid 
He flow cryostat under illumination with white light from a tungsten lamp. Resonant changes 
AI/I « 10~ 5 of the photocurrent I were detected using a current amplifier, lock-in detection 
and magnetic field modulation with 0.2 mT at 1.124 kHz. All spectra were obtained using 
a microwave power of 50 mW and are normalized to a microwave frequency of 9.749 GHz. 
Figure El shows the EDMR signal observed as a function of the applied magnetic field 
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for different angles between the [001] direction of the Si epilayer and the externally applied 
magnetic field Bq (rotation around [110]). To extract the spectroscopic information, the 
spectra were fitted with Lorentzian lines. In the range of B = 346.5 to 348 mT, anisotropic 

center at the Si/Si02 interface are observed with their characteris- 



tic anisotropy 



e . 


D 
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22j (dashed lines in Fig. |2J). Within the resolution of our experiment, 



these resonances remain unchanged in the strained layers. The dominant central line with 
g = 1.9994 in the unstrained sample becomes slightly anisotropic under strain and could be 
assigned to CB electrons [23j. In the unstrained Si layer (cf. Fig. 121(a)), also the charac- 
teristic two hf-split lines of P in Si with a separation of Ahf = 4.2 mT are easily resolved. 
The hf lines are isotropic as indicated by the black vertical lines with a center of gravity at 



g=1.9985, well known for unstrained c-Si (23, 



241- 



In contrast, the hf-split resonances become anisotropic for a strained epilayer as indicated 
by the black lines in Fig. 0(b) and (c). For Si on a Sio.93Geo.07 substrate (ej_=— 0.00199), 
they can be described with an isotropic hf interaction of 1.97 mT and an anisotropic g-factor 
which varies by A#=(1.46 ± 0.1) x 10" 3 from B Q || [110] to B || [001]. Figured (c) shows 
the results for the epilayer with e±=— 0.00729 obtained on a Sio.75Geo.25 substrate. Here, 
the isotropic hf splitting decreases to 0.94 mT, while Ag=(1.21 ± 0.06) x 10~ 3 . The small 
P-related hf splittings exclude significant segregation of P into the SiGe layers during the 
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growth process, which would lead to larger splittings 

For an isolated P donor atom in unstrained Si, the isotropic Fermi-contact hf interaction 
Ahf is given by 8/3/iB7r|-?/;(0)| 2 jll|, 0], where hb is Bohr's magneton and ^(O)] 2 is the 
probability amplitude of the unpaired electron wave function at the nucleus, giving rise to 
the doublet of hf lines separated by A^f = 4.2 mT in Fig. Ufa). To determine |?/>(0)| 2 , we 
first note that the cubic crystal field leads to the formation of a singlet ground state plus a 
doublet and a triplet of excited states instead of a six-fold degenerate ground state of the 
donor in Si. Only the fully symmetric singlet ground state has a non-vanishing probability 
amplitude at the nucleus, so that ip can be written as a superposition ip = ^^ =1 (l/v^6) < l ) i of 
the six valleys contributing to the donor. Here, each is a product of the corresponding CB 
Bloch wavefunction and a hydrogenic envelope-function. The probability of the unpaired 
electron at the nucleus (^(O)] 2 becomes 1/6| Y^=i ®i (0)1 2= 6 1^(0) | 2 , since due to symmetry 
$i(0)=$(0) for all i. Assuming that the only effect of strain is the change of relative 
population of the CB minima, we similarly find ^=E?=i( 1 /v / 2)^ and |^(0)| 2 =2|$(0)| 2 
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under high uniaxial strain, when only two CB minima contribute. Therefore, in the fully 
strained case, the hf interaction should be 1/3 of the unstrained case. In contrast, in Fig. El 
we already observe a reduction to 0.21 • A^-p, clearly below the 0.33 • Ahf limit obtained 



above. Based on group and linear elasticity theory, Wilson and Feher [11] have evaluated the 

H / 2 C \ 

analytical dependence of A H f(x) on the so-called valley strain x — t~ I 1 + „ ) ey, 

3A C V Ciil/I 



where E u = 8.6 eV 



28j is the uniaxial deformation potential, and 6A C = 2.16 meV 



is 



the energy splitting of the singlet and doublet state in the unstrained material. In Fig. |31 
a comparison of the prediction of Eq. (2) in Ref. \v\ (dashed line) with the hf splittings 
determined experimentally (full circles) clearly shows, that pure valley repopulation is not 
able to describe the experimental data for x > 0.07. An empirical treatment of additional 
radial redistribution effects as discussed in Ref. |30] would lead to 0.30 ■ A H f for X — ~ 89 
(x = 1), only a slight reduction of the repopulation limit and, thus, still at strong variance 
with the experimental data. 

An ab-initio calculation of hf interactions is necessary to clarify the situation, but is a 
demanding task. Already for the unstrained case, the delocalisation of the donor wave func- 
tion leads to problems in describing the magnetisation distribution correctly j^]. Whereas 
the energetic position of shallow levels as well as the geometries obtained by usual LSDA 
supercell calculations can be expected to be in accordance with experimental data, the hf 



interactions are generally overestimated 



12J, [32] . Recently, a Green's functions approach 
(LMTO-GF) has been shown to circumvent this problem of usual supercell approaches . 
By carefully analyzing the donor-induced resonance at the bottom of the CB, it becomes 
possible to describe the central-cell correction to effective mass theory by first principles 
with an accuracy that allows a prediction of super h ype r fine interactions for shallow donor 
states including the Kohn-Luttinger oscillations Q,!^. 

In the case of a strained host material, however, the situation becomes more complicated, 
since excited states are admixed to the former pure singlet ground state. An application of 
density functional theory (DFT) is only possible in combination with linear elasticity theory: 
Due to the applied strain, the symmetry of the P donor is reduced, and the resonance at the 
bottom of the CB is transforming according the a\ representation in the unstrained case, 
now shows admixtures of the b\ and b 2 representations of D 2 d symmetry. The location of 
the P donor atom in their nodal planes implies a correlation of these b\ and 62-like orbitals 
with the admixed doublet state. Since, furthermore, only one component of the diamagnetic 
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doublet state contributes to the singlet ground state under strain [ll|, it is reasonable to 
construct the spin-densities, which enter the self-consistent LSDA total energy calculations 
for a given valley strain, by n <T (r) = (l — <y(x))' n ai( r )~^ a (x) ' n w( r ); where a =|f)or \{). Here, 
a(x) is obtained from the strain-dependent admixture of the doublet states determined by 
linear elasticity theory (see Eq. (C6) in Ref. [11])- Figure El shows that the spin densities 
constructed this way allow a reasonable description of the pure valley repopulation effect 
since for an unrelaxed structure of an ideal Si crystal, the results obtained by Wilson and 
Feher jll| are nicely reproduced after the self-consistent cycle (cf. open squares in Fig. 

We are now able to take into account explicitely by first principles the strain of the 
Si lattice as well as the relaxation around the P donors within. For the optimization of 
the strained Si cells, we used the efficient self-consistent charge DFT-based tight binding 
(SCC-DFTB) approach Q- 

Optimisation of long slabs (up to 12 unit cells along the [001] 
direction) show an almost linear dependence of the compression along [001] as an answer to 
the tensile strain in the (OOl)-plane, effectively following linear elasticity theory as indicated 
in the inset of Fig. El This result confirms that linear elasticity theory remains valid in the 
complete regime investigated, even up to pure germanium as a substrate (x ~ -89). The hf 
parameters calculated with LMTO-GF under these assumptions (cf. open circles in Fig. EJ) 
already become smaller since the donor wavefunction becomes more delocalized as a result 
of the enhanced volume, exceeding the high-stress limit of 0.33 • Ahf obtained above. 

This tendency is enhanced, if local relaxation around the P donors is taken into ac- 
count: According to SCC-DFTB calculations on large, explicitly strained supercells with 
512 atoms, this relaxation is dominated by a slight reduction of the bond-length between 
the P donor and its nearest Si ligands by about 1%, nearly independent of the strength 
of the tensile strain in the plane of the Si epilayer and already present in the unstrained 
case. Re-calculating Ahf(x) f° r this geometry with the LMTO-GF code, we find a further 
reduction (full triangles in Fig. EJ). While the theoretically predicted absolute value Ahf is 
« 2 niT too large in the unstrained case without nearest neighbour relaxation, Ahf is now 
in good accordance with the experiment for all x- The hf interaction observed experimen- 
tally in the moderately strained P-doped Si layers can, thus, be explained by the increased 
volume of the unit cell together with a slight inward relaxation of the nearest Si neighbors. 
Since already such a small relaxation has a huge influence on the predicted relative hf split- 
tings for the strained material, the remaining discrepancy between experiment and theory 
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can easily be attributed to uncertainties due to the well-known flatness of the total energy 
surface in Si (35 1. According to our DFT calculation, the decrease of the central P-related 
hf interaction is accompagnied by a remarkable increase of the superhyperfine interaction 
with neighbouring 29 Si atoms. Comparative electron-nuclear double resonance (ENDOR) or 
electrically detected ENDOR (EDENDOR) [r| measurements could be used to verify these 
predictions. 

Finally, we briefly turn to the anisotropy of the g-f actor of the donor. The values observed 
here are somewhat larger than the values reported by Wilson and Feher |h| and those ob- 
served in strained two-dimensional electron gases J^fl . Note that a single- valley effect arising 
from an admixture of the doublet is not expected for uniaxial strain in 001] direction, dom- 



37 



381 ] . However, a 



inant in our samples. Rather, a pure repopulation effect is expected [n , 
detailed calculation of the ^-factor under strain including local relaxation appears warranted 
to understand the origin of the larger Ag observed in our samples. 

To summarize, we have presented an experimental and theoretical study of the hyperfine 
splitting of phosphorus donors in strained layers up to high strain levels of e±=— 0.00882. 
The splitting is reduced to values far below predictions published so far. Density functional 
theory demonstrates that repopulation of the CB minima by strain, the change of the unit 
cell volume, and a relaxation of the bond length between the P donors and the next nearest 
Si neighbours are required to account for the observed hyperfine interaction. Our results 
indicate that there exists no high-stress limit for the reduction of the P-related hf splitting. 
This work was supported by the Deutsche Forschungsgemeinschaft (SFB 631). 
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FIG. 1: (color online) Reciprocal space map of the 224 reflex of a 15 nm thick fully strained Si film 
on a Sio.7Geo.3 virtual substrate (log.-scaled iso-intensity plot, Cu Ka\ radiation). The straight 



line represents qu and q± predicted for relaxed Sii_ x Ge x with < x < 0.3 
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FIG. 2: (color online) Electrically detected magnetic resonance of P donors in relaxed and fully 
strained Si epilayers on three different Sii-^Ge^ virtual substrates as a function of the orientation 
of the [001] growth direction with respect to the magnetic field Bo for a rotation around [110]. To 
extract spectroscopic information, the spectra have been fitted using Lorentzian lines. The sum of 
these lines is shown in red. 
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FIG. 3: Hyperfine splittings observed for P in fully strained epilayers on Sii_ x Ge x substrates as 
a function x and the resulting valley strain. The black dots indicate the experimental data. The 



The DFT results for 



dashed line represents the behaviour expected by valley repopulation 
valley repopulation only are indicated by open squares, the results including the strain-induced 
change in the volume of the unit cell by open circles, and those additionally including relaxation of 
the P nearest neighbors by filled triangles. For the latter and the experimental data the absolute 
values ^4hf are shown, the remaining data refers to the relative scale. The insert shows the out- 
of-plane lattice constant a± of thin Si layers (thickness 12 unit cells) as a function of the in-plane 
lattice constant ay predicted by SCC-DFTB compared to linear elasticity theory indicated by the 
straight line. Pseudomorphic Si layers on pure Si and Ge substrates are shown by dashed lines. 
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